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ABSTRACT 


The optical path length from a satellite to the earth’s surface 
is s'^ongly dependent on the atmospheric pressure along the propa- 
gation path. Surface pressure can be measured by using a multicolor 
laser ranging system to observe the change with wavelength in the 
optical path length from the satellite to a ground target. The 
equations which relate surface pressure to the differential path 
lengths are derived and the accuracy of the pressure measurement 
is evaluated in terms of the ranging system parameters. The results 
indicate that pressure accuracies of a few millibars appear feasible. 
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I. INTRODUCTION 


Weather forecasting has benefited greatly from the rapid global 
coverage provided by meteorological satellites. Cloud cover, temperature 
profiles, pressure profiles and atmospheric wind speed are a few of the 
important meteorological parameters which can be measured from space. 

Cloud cover and temperature are readily observed using passive techniques, 
such as photography and radiometrv, while wind speed and pressure measure- 
ments require active lidar techniques.^ - Pressure is probably the most 

difficult of these parameters to measure accurately. Recently, Kalshoven 

2 

and Korb proposed a DIAL technique for the remote sensing of pressure 
and temperature. The system observes the pressure and temperature induced 
changes in the 7600 X oxygen absorption band. In this paper we describe 
another approach for measuring atmospheric pressure using a multicolor 
laser ranging system. The technique is based on the fact that the difference 
between the optical path lengths from a satellite to a ground target 
measured at two wavelengths is proportional to the atmospheric pressure 
at the target. 

The Spaceborne Geodynamic Ranging System (SGRS) currently under 

development at the Goddard Space Flight Center is a single color system 

3 

which operates at the 0.53 urn doubled YAG laser wavelength. The system, 
which will use retrof lector ground targets, is being designed to provide 
range resolution approaching one centimeter. Although accurate formulas 
have been developed to correct the range measurements for the effects of 
atmospheric refraction, an extensive weather station network is needed to 
provide the meteorological data required by the correction formulas. ” 
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To eliminate this problem, a multicolor version of SGRS using the fundamental 
(1.06 urn), doubled (0.53 urn) , and tripled (0.353 pm) YAG laser wavelengths is 
also being developed. ^ By processing the differential path length measure- 
ments using the approach described in the following section, this system 
could be used to measure surface pressure at the millibar level. In Section II 
we present the theory and derive the equations which relate surface pressure 
to the differential path length. In Section III the pressure measurement 
accuracy is examined. 
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II. THEORY 


The geometry of the satellite and target is illustrated in Figure 1. 

The optical path length is defined as the integral of the group refractive 
index along the ray path. Because the horizontal refractivity gradients 
are small, the one-way optical path length measured by a pulsed laser aystem 
is given by 
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n is the group index of refraction and 6 is given by Snell's law for a 
g 

spherically stratified medium (see Figure 1) 

nr cos 9 = n Q r 0 cos 9 Q (2) 


where n^ is the phase refractive index at the target. The atmospheric 

correction is the difference oetween the optical range R^ and the straight- 

line path length R . If n is expressed in terms of the group refractivity 
s g 


n = 1 + 10"S 
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Chen the atmospheric correction can be written in the form 
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The first tern is the velocity correction while the bracketed term is the 
difference between the geometric lengths of the ray and straight-line paths. 
The atmospheric correction can also be expressed as 

AC - SC + GC . (5) 
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SATELLITE 



Figure 1. Geometry of the satellite and target 


The spherical correction term SC corresponds to a spherically symmetric 

atmosphere while the gradient correction term GC includes the effects of 
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horizontal gradients. Marini and Murray have developed an accurate 
formula for calculating SC in terms of the surface pressure, temperature 
and water vapor pressure at the target during the satellite pass. The 
accuracy of their formula has been checked extensively by comparing it 
with ray tracing data. The results indicate that the error in Marini and 
Murray's formula is less than 5 mm for elevation angles above 20°. 

Recently, we developed an expression for GC in terms of the horizontal 
pressure and temperature gradients at the target.^ The accuracy of the 
gradient correction formula was also evaluated by comparing it with ray 
tracing data.^ because GC contributes less than 1 cm to the total atmos- 
pheric correction above 20° elevation, its effect will be neglected and 
we will restrict our attention to elevation angles above 20°. 

The atmospheric correction can be measured remotely from a satellite 
using either single-color or multicolor laser ranging systems. If the relative 
positions of the target on the earth's surface and the satellite are known, 
then the atmospheric correction can be calculated by subtracting the known 
distance between the satellite and target from the optical path length 
measured using a single-color ranging system. Because the straigh -line 
distance between the satellite and target must be known to within a few 
centimeters for this technique to be effective, it is probably not very 
practical. As an alternative, multicolor systems can be used to determine 
the atmospheric correction by calculating the difference between the optical 
path lengths measured at two wavelengths. 
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Although Marini and Murray’s formula can be used to accurately estimate 
the atmospheric correction at a single wavelength, some of the approximations 
made in deriving the formula are not valid for estimating the difference 
between the atmospheric correction at two wavelengths. Fortunately, it is 
easy to modify Marini and Murray's derivation to eliminate this problem. 

The approach is to evaluate the integrals in Equation (4) by solving Equation (2) 
for sin 0, substituting the result into Equation (4) and expanding the 
integrals in inverse powers of sin E. The resulting integrals are then 
evaluated by using the perfect gas law, law of partial pressures and the 
hydrostatic equation to obtain a suitable refractivity profile. The procedure 
is detailed in references 4 and 8. The most significant terms of the expansion 
for the round trip atmospheric correction difference are 
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and 


X » laser wavelength (urn) 
e ■ water vapor pressure at target (mb) 

P ■ surface pressure at target (mb) 

T * surface temperature at target. (K) 

9 ■ colatitude of target 

H ■ altitude of target above sea level (km) 

E ■ satellite elevation angle. 

In deriving Equation (6), we have neglected all terms contributing a few 

millimeters or less above 20° elevation to the total atmospheric correc' ion. 

The atmosr' eric pressure at the target can be determined by measuring 

Rp using a two-color laser ranging system and then solving Equation (6) 

for P. Because R. , is only weaklv dependent on water vapor pressure and 

temperature, accurate estimates of e and T are not required. This aspect 

is discussed in detail in the next section. Although the dominant variation 

2 

of Rp with respect to P is linear, the appearance of ? in the expression 
for B (Equation (8)) will introduce quadratic variations which cannot be 
neglected. Since the variation of K with respect to P is negligible, 
Equation (6) can be solved for P by using the quadratic formula 

P - -■» * ^ W /2 (U) 


where 
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III. ERROR ANALYSIS 


In a realistic system it will not be possible to determine the para- 
meters R^,,, E, e and T exactly. Therefore, it is importan.. -o assess the 
effect of errors in these parameters on the accuracy of the recovered 
pressure. An approximate expression for the root-mean-square (rms) 
pressure error is 
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c and are the rms errors in R.,, E, e and T, respectively. 
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The partial derivatives can be calculated by taking the derivatives of 
Equation (13) or by applying the chain rule to Equation (6) 
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The magnitude of 3P/3T is plotted versus elevation angle in Figure 2. 
Although the pressure sensitivity to temperature error is greatest at the 
lower elevation angles, its value above 20° is still quite small. At 20° 


elevation, a temperature error of 50 C would contribute less than one 
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Figure Pressure me isuremenc sensitivity to temperature errors 
as a function of the satellite elevation angle. 




millibar to the pressure error. Thus, only a crude estimate of the surface 
temperature is required. An accuracy of 20°C to 30°C is easy to obtain and 
should be adequate. 

The pressure sensitivity to errors in water vapor pressure is constant 

with respect to elevation angle (Equation (20)). A 10 mb error in water 

vapor pressure would contribute approximately one millibar to the surface 

pressure error. Since the water vapor pressure can approach 100 mb when 

the surface temperature and relative humidity are high, its effect cannot 

be ignored. It should be possible to predict the water vapor pressure to 

within 10 to 20 mb using a seasonally adjusted model supplemented with 

surface data obtained from existing weather stations. 

The elevation angle sensitivity (Equation (19)) is plotted versus E 

in Figure 3. The sensitivity is greatest at the lower elevation angles. 

At 20° a 1 mrad error in the elevation angle contributes almost 3 mb to the 

pressure error. The elevation angle can be determined from the laser 

pointing data. Since pointing accuracies of 100 yrad or less are easily 
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obtained with current technology, * the resulting pressure errors should 

be less than a few tenths of a millibar. 

The major error source for this technique is the differential path 

length measurement. The magnitude of the pressure errors caused by errors 

in the measurement of depends on the choice of wavelengths. The 

sensitivity 3 ?/ 3 R^ is plotted versus elevation angle in Figure 4 for the 

three possible combinations of the fundamental (1.06 4 m), doubled (0.53 ^m) 

and tripled (0.353 -m) YAG laser frequencies. The design of a space 

qualified multicolor ranging system based on these three wavelengths is 
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currently under study at the Goddard Space Flight Center. 
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Figure 3. Pressure measurement sensitivity to elevation angle 

errors as a function of the satellite elevation angle. 
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ELEVATION angle 




r 4-niltLvicy to differential path length 
Pre**ure MtelUte action an*i«. 


..fpiri as 



For the technique to be effective in measuring surface pressure at 
the millibar level, the differencial path length must be determined with 
an accuracy of at least a few millimeters. Streak tube detectors capable 


of providing the 5 psec time resolution required to measure path length 
differences at che millimeter level have been demonstrated in the 
laboratory. ^ However, the path length error is also dependent on the 
laser pulse width and shape, target dispersion effects, received signal 
strength and detection strategy. For example, if a maximum likelihood 
estimator is employed and the received pulse is a raised cosine 
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[O otherwise 

then the rms range error has the form 


J_ - CT/(2tt»'N) 
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(23) 


where C is the velocity of light, T i3 the full width of the optical pulse 
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and N is the average number of photoelectron3 per pulse. ’ “ Thus, it 
should be possible to measure the path length difference at the millimeter 
level if at lea t 100 photoelectrons are detected for each pulse and the 
received pulse widths are on the order of a few centimeters or less. 

The sensitivity to path length errors increases with increasing elevation 
angle while the sensitivity to elevation angle and temperature errors 
decreases with increasing elevation angle. The optimum elevation angle for 
minimizing the total pressure error can be calculated by substituting 
Equations (18) - (21) into Equation (17), and solving for the zero of 
iOp/5E. If the small temperature error effects are neglected, the optimum 


elevation angle is 


-ii/; 


E „ - sin _1 < 4.71 x 10" 3 (i(\ 1 ) - f(A,)) — — 
opt ) 12a 

R 12~* 


( 24 ) 


Because a_ will probably be on the order of lOOyrador less while o_ will 

*12 

be on the order of millimeters, the optimum elevation angle will normally 
be below 20°. The total pressure error for the 1.06 urn and 0.353 urn wave- 
length pair is plotted versus elevation angle in Figure 5 for several values 

of a . The elevation angle, water vapor pressure and temperacure errors 

*12 

were conservatively chosen to be 100 „rad, 20 mb and 20 C, respectively. 

Similar curves for the 0.53 um and 0.353 urn wavelength pairs are plotted in 
Figure 6. Based upon the results plotted in these figures, pressure accuracies 
of a few millibars appear feasible if the path length can be measured with 
an accuracy of a few millimters or less. 
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IV. DISCUSSION’ 


Surface pressure can be measured from a satellite by observing the 
change with frequency in the optical path length from the satellite to a 
ground target. 3ecause the recovered pressure is relatively insensitive 
to surface temperature and water vapor pressure, only rough estimates of 
these parameters are required. The major error source appears to be the 
differential path length measurement. To obtain pressure accuracies 
approaching a few millibars using the doubled and tripled YAG laser 
frequencies, the differential path length error must be no larger than one 
or two millimeters. Millimeter level accuracies appear feasible with 
current technology provided the received pulse widths are on the order of 
a few centimeters or less. If low dispersion retroref lector targets and 
centimeter laser pulses are employed, it should he possible to keen the 
received pulse widths within acceptable iimits. Adequate retroref lector 
targets could probably be manufactured and installed for less than one 
thousand dollars a piece, * but the approach is only practical for measurements 
over land masses. The major benefit of this pressure measurement technique 
will be realized it it can be adapced for use over the ocean. Although the 
sea surface is rough, at nadir it acts more like a mirror than a Lambertian 
reflector.'^ The local surface height of the ocean varies over a range of a 
few meters so that the diffuse components of the scattered oulses will be 
broadened to a few meters in length. Because of the large wavelength separation, 
the time-resolved speckle patterns in these diffuse pulses will be uncorrelated. 
However, there may be specular reflections which do not broaden the pulses 
and are highly correlated at the two wavelengths. If so. these signal 
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components could be used to determine the path length difference. There 

has been some work published on the reflection of CW and relatively long 

pulse (tens of centimeters and longer pulse widths) laser radiation from the 
14 

ocean surface, but to our knowledge, very little work has been done with 
centimeter pulses. This problem needs to be Investigated further before 
the utility of measuring surface pressure over the ocean using multicolor 
laser ranging systems can be adequately assessed. 
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